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ABSTRACT 

The Magnetic Field Satellite (Magsat) was launched on Oct. 30, 1979 into a nearly polar, sun- 
synchronous orbit, carrying a scalar magnetometer and a vector magnetometer. The satellite re- 
entered the earth’s atmosphere on June 11,1 980, having measured and transmitted more than th>xe 
complete sets of global magnetic field data. The data obtained from the mission will be used pri- 
marily to compute a currently accurate model of the earth’s main magnetic field, to update and 
refine world and regional magnetic charts, and to develop a global scalar and vector crustal magnetic 
anomaly map. 

This report describes the in-flight calibration procedure used for 39 vector magnetometer sys- 
tem parameters and gives results obtained from some data sets and the numerical studies designed 


to evaluate the results. 
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MAGSAT VECTOR MAGNETOMETER CALIBRATION USING 
MAGSAT GEOMAGNETIC FIELD MEASUREMENTS 


INTRODUCTION 

The Magnetic Field Satellite (Magsat) was launched from the Western Test Range at Vanden- 
berg Air Force Base in California into a nearly polar, sun-synchronous orbit on October 30, 1979. 

The orbit initially had a perigee of 352 km (219 miles) and an apogee of 578 km (359 miles). The 
satellite reentered the earth’s atmosphere on June 1 1, 1980, having measured and transmitted more 
than three complete sets of global magnetic field data. 

Magsat carried a scalar magnetometer and a vector magnetometer. The scalar magnetometer 
was an optically pumped, atomic resonance, dual-cell, self-oscillating cesium vapor magnetometer, 
capable of measuring the magnitude of the earth’s magnetic field with a resolution of about ±0.6 
gamma ( 1 gamma * 1 nanoTesla). The vector magnetometer consisted of three flux gate sensing ele- 
ments aligned along neatly orthogonal axes, capable of measuring each component of the earth’s 
magnetic field with a resolution of ±0.5 gamma. 

The data obtained from the mission will be used primarily to compute a currently accurate 
model of the earth’s main magnetic field, to update and refine world and regional magnetic charts, 
and to develop a global scalar and vector crustal magnetic anomaly map. Details concerning the 
mission and its objectives have been given by Langel, Reagan, and Murphy ( 1). 

Prior to the launch ot Magsat, the values of thirty -nine parameters associated with the vector 
magnetometer were determined in the laboratory at Goddard Space Right Center. It could not be 
assumed, however, that these values would remain unchanged after launch and under space condi- 
tions. Because of this a statistical procedure was developed for estimating the system parameteis 
from the data taken by the scalar and vector magnetometers in space. A computer program was 
written and the procedure was tested on measurements made by the magnetometers in the magnetic 
test facility at Goddard Space Right Center. The procedure has been and is continuing to be app'i *d 
to the real magnetic field data transmitted from Magsat. 
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This report describes the in-flight vector magnetometer calibration procedure and gives results 
obtained from some data sets and the numerical studies designed to evaluate the results. A FOR- 
TRAN listing of the calibration program can be obtained upon request. 

THE MAGSAT VECTOR MAGNETOMETER 

The triaxial fluxgate magnetometer flown aboard the Magsat spacecraft was capable of measur- 
ing the earth’s magnetic field vector components with a resolution of ±0.5 gamma. It was developed 
at NASA’s Goddard Space Flight Center by M H. Acuna and co-workers (2). A basic triaxial flux- 
gate magnetometer with ',-bit anaiog-to-digital converter had a dynamic range of ±2000 gammas. 
This range was increased to ±64,000 gammas in steps of approximately ± 1000 gammas by using three 
7-bit digital-to-analog converters as field offset generators, one for each sensor axis. Ground com- 
mand could select either of two redundant configurations (A or B) for the magnetometer electronics; 
an additional ground command could select one of two configurations for analog-to-digital converter, 
field offset generators, digital interface, and power converter. A feedback coil nulled the ambient 
field independently along each sensor axis. The reader is referred to a rep ort by Acuna et al (2) for 
details of the design of the Magsat vector magnetometer. 

The magnetic field values (y , , y,, y 3 ) associated with the three sensor axes are computed from 
the equations 

6 

y j * u j7 w j7 " ^ u jk w jk +w j8 (64 “ K j> 3+w j9 (64 " K j) 5+a j (f j“ b j):i* 1 -2-3 (1) 

k-1 

fj * (integer given by A/D converter) minus (zero level value) 

= fine measurement (-2048 < fj < 2047) (zero level value - 2048) 

Kj = integer from 7-bit D/A converter 
= coarse measurement (0 < Kj < 1 27) 
bj * bias in fj 
aj * scale factor 



Wj k = value associated with kth step of offset generator 
% 1000 X 2 k " 1 ; k 3 1 to 7 
u jk * 0 or * 

Uj 7 3 0 if the most significant bit of Kj is 1 
Uj 7 3 1 if the most significant bit of K ; is 0 
Uj k 3 the kth bit of Kj, Uj 0 being the least significant. 

For a detailed explanation of equation (1), the reader is referred to Acuna (3). Nominal values were 
determined for the 33 parameters Wj k , aj, bj (j 3 1 , 2, 3; k 3 1 . 2, . . . , 9) by various calibration pro- 
cedures prior to launch (2, 3). The values for the A-configuration are given in Table 1 . The cubic 
and quintic terms take into account small non-linearities in the vector magnetometer response 
function. 

An orthogonal reference frame was defined in the sensor assembly such that the angle between 
a magnetic sensor axis and the corresponding axis in the orthogonal frame was not more than 1 min- 
ute of arc. Let (x { , x 2 , x 3 ) be the components of the magnetic field vector in the orthogonal frame 
and Ojj the cosine of the angle between the ith sensor axis and the jth axis of the orthogonal system 
at the common origin of the two systems. The vector component; (x, , x 2 . x 3 ) in the orthogonal 
system and the magnetic field values (y, , y 2 , y 3 ) computed from ( 1 ) for the n on-orthogonal axes 
of the sensors are related by the equations 


y, 3 a u xj +sja 12 x 2 +a, 3 x 3 

(2) 

y 2 * Xj + ^ **23*3 

(3) 

y 3 3 a 3 ,Xi +s 3 a 32 x 2 +a 33 x 3 

(4) 

s, 3 sin Vj/vj 

(5) 

s 2 3 sin v 2 /v 2 

(6) 

s 3 3 sin v 3 /v 3 

(7) 

V 1 =a 12 x 2' r l 

(8) 
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v 2 * a 2l X 1 / r 2 <9) 

v 3 =a 32 x 2 /r 3 (,0) 

r l “**u + ®2* 1 12 X I + * 1 jj exp(5 1 h, 4 Xj > (11) 

r 2 * h^j ^ h 23 (6jh 2 4X.>) (1—) 

r 3 * h 3 j + ^2^32 x 3 + h 33 exp ( 62 ^ 34 X 3 ) (13) 

5 j * 1 if Xj > 0 (14) 

*0 if Xj •■* 0 
*-l if x, < 0 

5 2 »lifx 2 >0 (15) 

» 0 if x 2 ■ 0 
»-l if x 2 < 0 


aj, *(1 -a^ 2 -«j 3 ) ,/4 ( 16 ) 

^-d-^,-^)^ (17) 

033 »(1 -o§j -a§ 2 )* ( 18 ) 


Nominal values were determined in the laboratory at GSFC prior to launch for the 6 parameters a= k 
0 * k) and for the 1 2 parameters hj k . The values for the A-configuration are given in Table 1 . For 
further details see Acuna (2,3). 

The Sj , $2, and s 3 factors in equations (2), (3), and (4) and equations (5) through ( 1 5) by which 
they are computed are designed to model non-linearities in the vector magnetometer. 


PARAMETER ESTIMATION PROCEDURE 

The parameters to be estimated are w jk (j * 1. 2, 3; k * 1. 2 9), *■, bj (j * 1, 2. 3), a J2 . 

a 13* “21’ °23’ a 3i ’ **32' Xjj (i * 1, 2, . . . , n; j * 1. 2. 3) where n * the number of vector-scalar 

sets of measurements used in the estimation. If (x u , x i2 , x i3 ) are the true vector components at 
time Tj, then the corresponding true magnitude is given by 
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f i4 = (xfi + x i 2 2 +x? 3 ) ,/i . 


(19) 


If (f i( , f i2 , f i3 , f ;4 ) is a true set of values at time r ; , we let (f u . ? i2 , f i3 , f i4 ) be the corresponding set 
of field measurements at time Vj. If (y ;i , y i2 , y i3 ) are the true values of the magnetic field along the 
sensor axes at time t v then from equation ( 1 ) we have 


f ij = fry " u ij7 w j7 + £ u ijk w jk - w ij 8 (64 - K ii> 3 


k-1 


“ w i,9 ^ ” K ij) 5 • / a j + bjii * 1 to n;j 3 1, 2, 3; (20) 

where the y^’s are computed from equations (2) through ( 1 8). We abbreviate equations [( 20), (2) 
through (18)1 by 


and equation ( 1 9) by 


fjj*«j(Xj,p) ; j 3 1, 2, 3; 
f i4 * g 4 (x t ). 


( 21 ) 


( 22 ) 


where gj 0 * 1 . 2, 3, 4) are function symbols, x i 3 (x n , x i2 , x i3 )’ , using' as a symbol for matrix 
transpose, and p is a column vector with elements Wj k (j 3 l,2,3;k 3 1 to 9), aj, a 2 , a 3 , b,, b 2 , b 3 , 
a 12’ a i3» °2i • CL i3' a 31 ’ ®32 ™ *** e order given. We further abbreviate for i 3 1 to n 


f i“( f it’ f i2’ f i3’ f i4^’ 

(23) 


(24) 

f t “ t^ii , f i2 , f i3 , . 

(25) 

A A A A A 

f i2’ f i3’ f i4> • 

(26) 

f ij 3 gj(S? i , ?),j 3 l,2,3. 

(27) 

f i4 =g 4 (jl i ). 

(28) 

fy * gj($i, p),j= 1,2,3, 

(29) 

f i4 - g 4 (X|). 

(30) 
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where x . and p are guesses and Xj and p are estimates for the true values of the parameters x s and p. 
By an estimate of a parameter we mean a value obtained with a statistical estimation procedure. 

We assume that all sets of vector-scalar measurements f. have the same error distribution with 
0 means, 0 correlations, and standard deviations (Oj , o,, o 3 , o 4 ), and that the a priori values of p k 
(k * 1 to 39) are samples of size 1 with error distributions whose means are all 0 and whose standard 
deviations are (k a 1 to 39). 

Suppose the true value associated with the ith measurement is <I/ jt that is given as a function 
of the complete set of true parameters 9 * (0 1 , 0 2 , . . ■ 9 m ) by i// i 3 g^O), and that errors in measure- 
ments of \l > j have mean 0 and standard deviation o ; . Then the “weighted sensitivity matrix” T is de- 
fined as the matrix whose element in the ith row and jth column is (d 'bJb9 j )lo l . Since we will have 
only a guess for 9 , ^ * g^ 9 ) and T is approximated by T. We see that T for our estimation will have 
4n rows and 3n + 39 columns. 

Since many of the partial derivatives for our calibration problem are 0, it will be convenient to 
partition I\ 

r-lRSl. (31) 



Rj 0 ... o" 


Si 

R» 

■ • O 

■ -J" 

• o 

, s- 

s : 


b 0 ... R n _ 


_S„_ 


and O is a 4 x 3 matrix all of whose elements are 0. R; S I r ijkl (i • 1 to n) are 4 x 3 matrices with 
r ijk 0’ * 1 to 4 f k * 1 to 3) as the element in the jth row and kth column with 

r ijk * )iGy 

TheO’sare zero submatrices of R because dfjj/dx hlc - 0 unless i * h. S { = [s i j k ] (i * 1 to n) are 4x 39 
matrices with Sy k (j = 1 to 4, k = 1 to 39) as the element in the jth row and kth column with 

SUk-Ofy/aPkWoj. (33) 

Note that all the elements in the 4th row are 0 for i = 1 to n. 
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Let z be a column vector of the true values of the estimated parameters and 12 a diagonal matrix 
whose diagonal elements are the variances for the elements of z, the a priori approximation of z. We 


define 


V 


V 

d : 




v d- - 


* 

» i 

d i3 

<*« 




Thus d is the weighted vector of residuals. Also 


< f y 


f ii ) /<y j • 


dy-ffy-FiP/Oj. 

The weighted least squares estimate of z is 

z-2 +K" 1 (?' d +I2" 1 (z -z)l, 

K*r'r+n- 1 . 


(34) 


(35) 


(36) 

(37) 


For the Magsat calibration problem, we assume infinite standard deviations for a priori values 
of Xy . This means the first 3n diagonal elements of I2” 1 are 0 and the remaining 39 are 1 /cj k 
(k * 1 to 39). Thus J2” 1 (z - z) is a column vector whose first 3n elements are 0 with the remaining 
elements given by 

t k * (p k - P k )/« k : k = 1 to 39. (38) 


We note that 


t k = (p k - P k )/<*J k 


and define a column vector t with 39 elements whose kth element is given by (39). We also define a 
39 X 39 diagonal matrix J with kth diagonal element equal to 1 /w£ . Thus 


n- 1 = 




(40) 


where Oj and 0 ; are 3n X 3n and 3n X 39 zero matrices. For our calibration problem we also have 
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(41) 


z * 



fC can be written in partitioned form as 



R' 



K* 


(R SI + 



S'_ 


O J 



R'R 

R'S 


A 

B 

* 
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S' R 

S'S + J_ 


_B' 

C _ 


where A*R'R,B a R’S, and C ~ S' S + J. Note that A, C, and K are symmetric. We also have 



y 1 

o 

o 

1 


V 

A* 

0 A 2 ... 0 

, B* 

*2 

“ ! 

* 1 


0 6 ... A„_ 


."■J 


Aj * R' R if Bj * R' Sj (1*1 to n) 


C* £ S'S t + J 

i-1 


(43) 

(44) 

(45) 


Aj is 3 x 7 . Bj is 3 x 39, C is 39 x 39, A is 3n x 3n, B is 3n X 39, K is (3n + 39) x (3n + 39). 
Since n is typically several thousand, a direct inversion of K is impractical. The inversion problem 
is manageable, however, if we take advantage of the structure of A in a partitioned fonr of the 
inverse. 

We can now write the bracketed part of equation (36) in the form 


r'd + sr 1 (z-z)* 


R 1 


“ - 


“ — 

d + 

O 

3 

0 

v m 

t 


_7_ 


(46) 


8 


with the definitions 


0=R'd. (47) 

7 *S'd+t, (48) 


where 0 is a column vector with 3n components, 7 is a column vector with 39 components, and t is 
defined by equation (39). 

Equation (36) becomes 


A 

X 


X 


A B 


V 


a 


+ 




A 

_P_ 




B' C 
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and we recoil that A’ * A and C' * C. 

If we express the inverse in f he form 


A 

B _ -> [*D 

e" 

JB' 

cj |je' 

G_ 


(49) 


( 50 ) 


and recall the well-known formulas for the inverse of a partitioned matrix, we have 

H * -A~ l B, 

G * (C + B' HT 1 , 

E * FG. 

D* A " 1 + EH'. 

Partitioning H, E, and D as 



* 


V 



D i: 

l 

C 

a 

H* 

« 2 

, E * 

E 2 

, D = 

■>21 

D 22 • 

■ D, n 




_ E n_ 



D „2 • 

O 

3 

\ 


(51) 

( 52 ) 

(53) 

(54) 


( 55 ) 
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D[* - D- f , i = l to n t j " I to n, and using (43) we have 

= -A " 1 Bp i - 1 to n. 
n 

G * (C + y B; Hj)- 1 . 

l-l 

Ej = H ; G; i = 1 to n, 

Dj| * Aj -1 + Ej H' ; i 3 1 to n, 

Djj = E i H- : i = I to n - 1 : j * i + 1 to n. 

Equation (49) can be written as two equations 

x » x + D 0 + E 7 , 

p * p ♦ E' 3 + G 7 . 

From equation (47) we obtain 



V 


^il 

0= 

h 

. 

fa 


0 n 


A 


L. -i 

1 



We can now write (61 ) and (62) in the forms 


< 56 > 

(57) 

(58) 

(59) 

(60) 


(61) 

(62) 


(63) 


x i *x i + E i7 + £ Dyfy i-i ton 
j-1 


n 

p-p+G7+ J e; 3 j. 

i«l 


(64) 


(65) 


We note that E t is 3 X 39, Dy is 3 X 3, G is 39 X 39, G' = G, and recall that D' } = Dy for all i’s 
and j’s. 
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The covariance matrices associated with the estimates x i (i = 1 to n) and p are given by 


cov(x., Xj)* D ir 

(66) 

covfXj, p)= E|, 

(67) 

cov(p, p) ■ G. 

(68) 


Correlation matrices can be obtained from the covariance matrices. 

RESULTS OF A CALIBRATION 

The calibration results given below were obtained by exercising the computer program on a 
small subset of the data recorded on Nov. 5, 1979. The number of vector-scalar sets used in the cali- 
bration was 1 247, approximately equally spaced in time. Table 1 gives the ground-based calibration 
values for the system parameters along with the standard deviations used to compute the weights for 
these a priori values in the least-squares algorithm. Except for the Oy parameters, these standard devi- 
ations are larger than the realistic ones for the ground-based calibrations. This allows for the fact that 
the true values of the parameters may be significantly different after launch than before launch. The 
standard deviations for the a pnori values of a 12 , a i3 ,a 21 ’ ^3’ *33 are held at the ground-based 

calibration values because of strong correlations between the errors in the least squares solutions for 
these parameters. Table 1 also shows the differences between the ground-based calibration values 
and the values estimated from the Nov. Sth (79) data. The final column of Table 1 gives the standard 
deviations of the estimated parameters. The three parameters w 25 , w 26 , and w 27 enter the least 
squares procedure only as the linear combination w 27 - w 26 - w 2J , due to the fact that the magnetic 
field components along the second sensor axis had maximum and minimum values of 10228.5 and 
-1 5246. 1 . To obtain a solution, w 25 and w 26 were held fixed by giving them very small a priori 
standard deviations. 

Table 2 gives 5 examples from the 1247 sets of vector-scalar measurements used in the least- 
squares calibration for Nov. 5, 1979. Kj and (i * 1 , 2, 3) are the coarse and fine values recorded 
for the three magnetic field components and f 4 is the field magnitude measurement. The Kj’s are 
error-free and are not considered as measurements in the least-squares estimation. procedure. 
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Table 1 


Parameter 

Pte-launch 

Value 

Standard 

Deviation 

W 11 

999.7 

20 

W 1 2 

1998.7 

20 

W 13 

3996.7 

20 

w 14 

7993.3 

20 

W 15 

15986.1 

20 

w 16 

31973.8 

20 

w 17 

63947.5 

20 

w 18 

.963 E-5 

.32 E-5 

* 19 

-.603 E-9 

.32 E-9 

W 21 

1004.3 

20 

w 22 

2008.7 

20 

w 23 

4016.8 

20 

w 24 

8037.6 

20 

w 25 

16069.6 

1 E-l 2 

w 26 

32143.9 

1 E-l 2 

w 27 

64287.9 

20 

w 28 

.436 E-4 

J2E4 

w 29 

-.649 E-9 

32 E-9 

W 31 

997.3 

20 

w 32 

1995.1 

20 

w 33 

3989.8 

20 

w 34 

7979.6 

20 

w 35 

15958.5 

20 

w 36 

31918.2 

20 

w 37 

63844.2 

20 

w 38 

-.508 E-5 

.32 E-6 

w 39 

.167 E-8 

.32 E-9 

bt 

4.0 

20 

b 2 

13.0 

20 

b 3 

12.0 

20 

a l 

.98248 

1 E^ 

a 2 

.99226 

1 E4 

a- 

.98322 

1 E4 

La l2 

-228.5 

6.5 

i3 

168.0 

2.1 

*- a 21 

-612.3 

6.5 

L *23 

443.3 

2.1 

^ a 3l 

483.1 

2.1 

L °32 

589.9 

6.5 


Estimated 

Correction 



-0.4 
- 1.1 
-1.9 
-4.5 
- 10.2 
-0.4 E-6 
-0.5 E-l 1 
0.7 
-0.3 
-0.5 
-2.9 
0.0 
0.0 
0.1 

-0.7 E-5 
-0.3 E-l 2 
0.5 

1.0 

1.9 

3.0 

4.8 

8.8 

18.5 

0.6 E-5 
0.3 E-9 
6.8 

4.6 

0.1 

0.4 E-4 
-2.1 E-4 
-1.7 E-4 
-16.8 

1.0 

8.0 

2.7 

1.0 

28.0 


Standard 

Deviation 


.05 

.05 

.06 

.08 

.14 

.29 

.25 

.31 E-5 
.32 E-9 
.16 
.17 
.25 
.34 
1 E-l 2 
1 E-l 2 
.46 

.31 E-4 
.32 E-9 
.03 
.03 
.04 
.05 
.09 
.14 
.12 

.31 E-6 
.26 E-9 
.24 
.48 
.19 
.5 E-4 
.9 E-4 
.4 E-4 
4.2 

1.5 

4.7 

2.0 

1.5 
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Table l 

(unestimated parameters) 


h ll 

-35.28 

h l2 

1.46 E-4 

h 13 

-3.78172 

h t4 

-4.61641 E-5 

h 21 

-116.45 

h 22 

3.35 E-4 

h 23 

-10.832 

h 24 

-3.33 E-5 

h 31 

-127.1 

*32 

3.2 E-4 

*33 

-13.3367 

h 34 

-2.9167 E-5 


40 y * the complement in arc seconds of 
. the angle associated with the direc- 
tion cosine ay 


Table 2 


Kj 

fl 

mm 

h 

k 3 

h 

^4 

55 

-900 

62 

-655 

14 

-563 

51453.9 

32 

- 376 

59 

-417 

65 

522 

32740.9 

74 

855 

49 

812 

99 

-462 

38780.4 

45 

754 

61 

613 

71 

719 

19981.9 

54 

136 

73 

-1012 

106 

968 

44726.4 


13 




















In the remainder of this section w will present some computer results obtained for these 5 


vector-scalar sets as well as for all 1 247 sets. 

Table 3 gives the distribution of the absolute values of 64-K j for the 1 247 data set . 

Table 3 


Number 

64-K, 

64-K j 

64-K 3 

-0 

21 

46 

11 

- 1 

40 

177 

14 

>140 

68 

287 

22 

>340 

132 

437 

65 

>7405 

367 

300 

131 

>15401 

603 

0 

339 

>314<63 

16 

0 

665 


The values iy t , y 2 , y 3 ) of the magnetic field along die three sensor axes are obtained by sub- 
stituting the pre-launch values of 33 parameters and the data in Table 2 into equation (1 ). The re- 
sults are given in the first 3 columns of Table 4. The fourth column of Table 4 gives f 4 - f c , where 
f c is the approximate vector magnitude computed from the components y, , y 2 , y 3 , using a non- 
Pythagorean formula. 


Table 4 


y\ 

y 2 

>3 

f 4 -f c 

-9881.2 

-2671.5 

-50437.0 

15.1 

-32347.4 

-5447.8 

1507.1 

1.5 

10828.3 

- 14274.7 

34452.8 

2.5 

-18247.6 

-2417.6 

7685.7 

4.0 

-9862.3 

8027.5 

42840.6 

10.0 
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When the values of y ,, y 2 , and y 3 from Table 4 and the prepunch values of a,-,, a, 3 , a,, , a, 3 , 
Oj j, and a 32 are substituted into equations (2) through (18), the vector components (.;, , x,, x 3 ) in 
the orthogonal system can be computed. They are given in the first three columns of Table 5. The 
fourth column of Table 5 gives f 4 - f , where f is the vector magnitude computed from the compo- 
nents x j , x 2 , x 3 , using the Pythagorean formula ( 1 9). Differences between the fourtn columns of 
Table 4 and Table S are due to the fact that Sj , s,. and s 3 differ slightiy from 1 . 


Table 5 


*1 

*2 

*3 

f 4 - f 4 J 

-9843.0 

-2592.1 

-50421.7 

i 

15.1 j 

-32354.8 

-5537.8 

1567.1 

3.1 

10784.8 

-14316.9 

34387.2 

2.0 

- 13256.7 

-2486.7 

7721.4 

4.2 

-9888.6 

7906.3 

42886.6 

10.1 


The values of x 3 , x 2 , and x 3 in Table 5 are used as first guesses in the least squares estimation 
procedure. 

Table 6 gives the values of x j - Xj , x, - x 2 , and x 3 - x 3 in the first three columns, where 
(Xj , x 2 , x 3 ) are estimated values and (Xj , x 2 , x 3 ) are the initial guesses in Table 5. Column 4 of 
Table 6 gives the values of f 4 - f 4 , where f 4 was computed from x,, x 3 ) using the Pythagorean 
equation. We note that the estimated field vectors have magnitudes much closer to the measured 
magnitudes than the guesses (Xj , x,, x 3 ) computed from the vector field measurements , fU, fg) 
and the a priori parameter values. The means and root-mean-squares of the 1 247 values of x t - x , . 
x 2 - x 2 , x 3 - x 3 , and f 4 - f 4 are given ; n Table 7. 

Let (Xj , x 2 , x 3 ) be the values computed by substituting the estimated 39 parameters in equa- 
tion (1) through (18) and solving for (X[ , x 2 , x 3 ), and let f 4 * (xj + \\ + x^)^. Table 8 shows the 
differences between these estimates for (x t , x 2 , x 3 , f 4 ) and the least squares estimates (x, , x 2 , x 3 , 
f 4 ). 
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Table 6 


x l -X t 


x 3- x 3 

ns 

5.1 

5.1 

14.6 

-0.5 

2.0 

5.6 

2.5 

0.3 

8.9 

0.6 

-5.9 

-1.0 

3.7 

5.6 

- 0.7 

-0.2 

5.6 

4.1 

-11.0 

-1.0 


Table 7 



“m A 

X, -X, 

KS9 

ESS 

ns 

MEAN 

7.0 

3.2 

1.7 

-0.001 

RMS 

3.0 

1.8 

9.8 

0.52 


Table 8 


A V 

x,-x, 

ESS 

mmm 

KBS 

0.1 

0.0 

0.4 

-0.4 

-0.2 

-0.1 

0.0 

0.2 

-0.2 

0.3 

-0.6 

-0.7 

0.1 

0.0 

0.0 

-0.1 

0.1 

-0.1 

-0.7 

-0.7 


The 1 247 values of f 4 - f 4 have a mean of -.002 and an RMS of .90, which compare very well 
with the values of -.001 and .52 in Table 7 for f 4 - f 4 . Thus least squares estimation is not needed 
for every day’s data if a loss of about 0.5 gamma in accuracy can be accepted. We must consider 
that on the GSFC 360/91 , full estimation for 1 247 data sets required 8 minutes, while computation 
of (Xj , x 2 , x 3 ) from equations (1) through (18) required 0.4 of a minute, a ratio of 20 to 1. 

The first three columns of Table 9 give the residuals fj - fj (j * 1 , 2, 3 \ where the f^’s are the 

A 

fine measurements given in Table 2 and the fj are obtained by substituting the estimated parameters 
in equation (20) and computing (y 2 , y 2 , y 3 ) for (20) by substituting (Xj , x 2 , x 3 ), a 12 , a 13 , cL,j , 
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a-j j , and a 32 in equations (2) through (18). The fourth column of Table 9 repeats the fourth 
column of Table 6. The 1 247 value sets of the numbers given for a sample of 5 in "U 9 have 
means and RMS’s given in Table 10. Table 10 is an indication of the goodness-of-fit of the solution 
to the measurements. 


Table 9 


w A 

■JBM 

• A 

f 3 -f 3 

mm 

-0.1 

-0.0 

l 

-0.5 

0.2 

0.0 


0.3 

0.2 

-0.3 


-1.0 

-0.1 

1 

mm 

-0.2 

-0.2 

0.1 

■91 

-1.0 


Table 10 


_ A 

f 2 -f 2 

w A 

f 3 -f 3 

f 4 -f 4 

-.001 

-.0009 

.0002 

-.001 

0.19 

0.06 

3.32 

0.52 


FURTHER CALIBRATION RESULTS 

Knowing that under space conditions the calibration parameters determined pre-launch would 
vary, it was necessary to determine criteria for a first in-flight calibration, to be used as a control, to 
see any trends and shifts in subsequent calibrations typifying changes in the instrument. November 5th 
data was chosen for the control calibration under the criteria that ( 1 ) the instrument should have had 
time to shift mechanically, due to stresses in a space vacuum system, into some “permanent” position 
after three days in orbit, and (2) data were taken continuously throughout the day. producing a glob- 
ally distributed data set for sampling. 

The November 5th calibration (see Table 1 1) shows changes from pre-launch calibration (espe- 
cially in the z-axis), which illustrate the changes in the instrument. Calibrations were then made at 
four day intervals beginning November 2. 1979 and ending January 21,1 980 (excluding January 1 , 
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Table 1 1 


Parameter 

Pre- Launch 
Value 

November 5th 
Value 

Estimated 

Correction 

Standard 

Deviation 

w u 

999.7 

999.6 

-0.1 

0.14 

"12 

1998.7 

1998.3 

0.4 

0.14 

w 13 

3996.7 

3997.0 

0.3 

0.15 

w 14 

7993.3 

7993.2 

-0.1 

0.16 

w 15 

15986.1 

15985.9 

-0.2 

0.24 

w 16 

31973.8 

31973.4 

-0.4 

0.43 

w 17 

63947.5 

63940-0 

- 7.4 

0.39 

"18 

.963 E-5 

.882 E-5 

-.803 E-6 

.98 E-6 

W 19 

-.603 E-9 

-.611 E-9 

-.829 E-ll 

.10 E-9 

W 21 

1004.3 

1005.3 

1.1 

0.42 

"22 

2008.7 

2009.0 

0.4 

0.41 

w 23 

4016.8 

4017.4 

0.5 

0.43 

w 24 

8037.6 

8036.7 

-0.9 

0.43 

w 25 

16069.6 

16068.8 

-0.8 

0.88 

"26 

32143.9 

32143.2 

-0.8 

0.88 

"27 

64287.9 

64288.7 

0.8 

0.88 

W 2 B 

.436 E-4 

.436 E-4 

-.910 E-8 

.10 E-5 

"29 

-.649 E-9 

-.649 E-9 

-.247 E-13 

.10 E-9 

"31 

997.3 

997.7 


0.10 

"32 

1995.1 

1996.0 


0.10 

"33 

3989.8 

3991.6 


0.11 

"34 

7979.6 

7982.7 


0.12 

"35 

15958.5 

15963.1 

4.6 

0.16 

"36 

31918.2 

31926.7 

8.5 

0.22 

"37 

63844.2 

63862.6 

18.4 

0.24 

"38 

-.508 E-5 

-.217 E-5 

0.291 E-5 

.96 E-6 

"39 

.167 E-8 

.173 E-8 

0.628 E-10 

.10 E-9 


4.0 

7.8 

3.8 

0.35 

b 2 

13.0 

13.9 

0.9 

0.44 

b 3 

12.0 

12.7 

0.7 

0.31 

a l 

.98248 

.98248 

0.0 

.50 E-5 

a 2 

.99226 

.99226 

0.0 

.50 E-5 

a 3 

.98322 

.98322 


.50 E-5 

£*12 

-228.5 

-228.3 

-0.2 

0.99 

£*13 

168.0 

168.0 


0.033 

^•*21 

-612.3 

-612.2 

-0.1 

1.0 

£■*23 

443.3 

443.3 


0.033 

£*31 

483.1 

483.1 

0.01 

0.033 

£*32 

589.9 

595.6 

5.7 

0.95 
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1 980, at this writing). The first 29 plots following this discussion show graphically the changes from 
the November 5th calibration for those parameters which changed more than 0.1 gamma (or 0.1 arc* 
seconds, according to which parameter is being considered). A key to the symbols used in the plots 
can be located on the page before the plots. Two sets of parameters, w I6 - w J7 , and w 2J - w, 6 -w, 7 
are plotted together since correlations between these parameters are so large that the parameters are 
believed to be mathematically inseparable. 

Examination of the calibration parameter plots shows trend-line changes with time, usually in 
one direction (i.e., non-oscillating), with the exception of several days calibrations. Upon further 
examination, these calibration exception* are usually found common in several plots. One reason 
for the offset values is the status of die scalar sensors, A and B. It is found that a different calibra- 
tion will result for each of the following cases: ( 1 ) sensor A on only, (2) sensor B on only, or (3) 
both sensor A and B on. Examination of the plots implies this to be the case in most instances. For 
those instances where the sensor status does not apply, other reasons such as diversion in the data 
distribution, or temperature variation of the sensors and platform may be part of an explanation, 
but no conclusions have been formalized at this time. 

The second set of plots describe some statistics about die relationship between the scalar values 
and their associated vector magnitudes, for samples of data at four day intervals about the date of 
the calibration. The data samples had the calibration in question applied to them, and the mean, 
standard deviation about the mean, and the maximum of the differences between scalars and vector 
magnitudes were determined. Study of these plots shows conclusively the dependence of the calibra- 
tion and the calibrated data on the scalar sensor status. 
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Key to Symbols used in Plots 


A: Scalar Sensor A on only 

B: Scalar Sensor B on only 

C: Scalar Sensors A and B both on 

M: Mean of difference values between sedars and 
vector magnitudes 

R: Standard deviation about the mean of difference 

values between scalars and vector magnitudes 
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DEVIATION OF (W17-WI6) ESTIMATE FROM NOVEMBER 5TH ESTIMATE 
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DEVIATION OF W25 ESTIMATE FROM NOVEMBER 5TH ESTIMATE 
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DEVIATION OF W27 ESTIMATE FROM NOVEMBER 5TH ESTIMATE 
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DEVIATION OF (W27-W26-W23) ESTIMATE FROM NOVEMBER 5TH ESTIMATE 
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SUMMARY 


In summary, the results we have described indicate first that the calibration parameters of the 
MAGSAT vector magnetometer, particularly those associated with the third (z) axis, changed signifi- 
cantly between the pie-launch calibration and the first post-launch calibration. This can possibly be 
attributed to slight mechanical shifts due to the difference in stresses between a vacuum and non- 
vacuum environment. Second, some of the calibration parameters underwent small and slow system- 
atic changes during the mission lifetime; and third, the calibration parameters determined from inter- 
comparison with die scalar instrument are slightly different (*• 1 tamma) between the cases when both 
sensor A and sensor B of the scalar magnetometer are operating, when only sensor A is operating, and 
when only sensor B is operating. There are three possible causes for differences in vector calibration 
with respect to the scalar sensor configuration: 

1. The change in null-zone configuration of the scalar instrument results in a bias in die distribu- 
tion of vector directions visible in the calibration. This could cause imprecise determination of 
some calibration parameters, particularly those defining cross-coupling between axes of the vec- 
tor instrument. 

2. With only one sensor operable, heading errors on the scalar magnetometer are increased. Such 
errors will be geographically preferential, which means they will also occur at systematic direc- 
tions of the field with respect to the vector magnetometer. That is, for certain field directions 
the scalar magnetometer would give an erroneous reading of up to 1 .5 gammas, with a corre- 
sponding change in the determination of the vector magnetometer calibration parameters. 

3. An interference of the scalar magnetometer with the vector magnetometer which changes with 
die scalar magnetometer configuration. Such an interference was experienced during pre-launch 
magnetic testing. In that case, the R. F. excitation signal from the scalar instrument caused a 
sensitivity shift in the vector instrument. This was most pronounced in the second (y) axis, 
still noticeable in the third (z) axis, and least discemable in the first (x) axis. Subsequent radio 
frequency shielding reduced the effect to less than one gamma at full scale. 
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In addition to knowing that the calibration parameters differ with the scalar configuration, we 
know that the calibration constants determined under a specific scalar configuration result in a larger 
mean difference when applied to data with a different scalar configuration. This would seem to rule 
out the first hypothesized cause for the difference. We need, however, to distinguish between causes 
2 and 3 because if cause 2 is correct then we should utilize the more accurate calibrations fi _>m times 
when both sensors A and B were operating. If, on the other hand, radio frequency interference dif- 
ferences are affecting the measurements, then the calibration actually does change and the calibration 
parameters should be those derived from the scalar data with the same configuration that was in effect 
at the time of the data being calibrated. 

Examination of the calibration parameter plots indicates that all axes are affected by scalar con- 
figuration changes. This seems to indicate that cause 2 is operable. However, more analysis is required 
before a final verdict is in. 
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